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A theory for electronic Raman scattering in the cuprate superconductors is presented 

with a specific emphasis on the polarization dependence of the spectra which can infer the 

symmetry of the energy gap. Signatures of the effects of disorder on the low frequency and 

' ^ ■ low temperature behavior of the Raman spectra for different symmetry channels provide 

C—, ' detailcd information about the magnitude and the phase of the energy gap. Properties of 

Q ' the theory for finite T are discussed and compared to recent data concerning the doping 

O dependence of the Raman spectra in cuprate superconductors, and remaining questions 

are addressed. 
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1. Introduction 
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_ ^ Interest in electronic Raman scattering has grown considerably since the first the- 

oretical analysiad in 1961 and the subsequent observation of the Raman effect in 

C^- [ superconductors seventeen years ago.H During the early- to mid- eighties the amount 

of available data was limited to a few A- 15 superconductors while the theory was 
elarified in many aspects.B With the discovery of the high temperature supercon- 
ductors, this situation dramatically changed as the amount of data grew on these 
systems. The theory of the Raman effect in superconductors was completed by 1991 
for both clean and disordered s-wave superconductors,Q but this theory could not 
capture many features shown in experiments on the cuprates. At present, experi- 
ments on the Raman effect exist for nearly all "high" temperature superconductors: 
materials which are electron or hole-doped and materials with different number of 
CuC>2 planesB 

Recently attention has turned towards uncojiventional superconductivity can- 
didates to describe the pairing in the cuprates.H Subsequently, due to the strong 
symmetry dependence of the observed spectra, i.e., the characteristic features of 
light scattering for different incident and scattered polarization orientations, in all 
high T c compounds, electronic Raman scattering in unconventional superconduc- 
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tors has grown to be of considerable interest in light of idcntifying the symmetry 
of the energy gap in high temperature superconductors. Indeed, the amount of at- 
tention lavished on this area has been remarkable and has provided a large amount 
of detailed information towards understanding the mechanism of pairing in these 
materials. 

The importance of Raman scattering can be related to its ability to probe ex- 
citational dynamics on regions of the Fermi surface rather than being restricted 
to measuring averages over the Fermi surface. The symmetry dependence of the 
spectra has augmented an understanding of the magnitude and symmetry of the 
energy gap, and makes it as effective a probe as photo-emission has proven to be 
to determine A(k). Moreover Raman scattering is a bulk probe of a material due 
to the long wavelength of the exciting laser light, does not suffer appreciably from 
surface effects, and has extremely sharp energy resolution. 

Simple considerations of the transformation properties of the light scattering 
amplitude were used in Ref. 7 to demonstrate how the light polarization orienta- 
tions selectively probe excitational dynamics of regions of the Fermi surface or the 
Brillouin zone, and subsequent work has clarified this picture considerablyE'uEj In 
the normal state, the light scattering cross section provides information concerning 
the scattering rate of electrons in certain k-space regions. B\ g scattering geometry 
[{x — y)(x + y) orientation of incoming (scattered) light polarizations, respectively] 
probes excitations along the Brillouin zone axes while Bi g \(x)(y)\ probes the diag- 
onals. Thereby k-dependent scattering processes are measured simply by rotating 
the polarization orientations of the incoming and scattered photons. In the su- 
perconducting state, a k dependent energy gap can be mapped out. The direct 
coupling of the Raman vertex to an anisotropic energy gap A(k) leads to a strong 
polarization dependence of the Raman spectra in the superconducting state. For 
an energy gap of d x 2_ y 2 symmetry, which is minimal along the zone diagonals and 
maximal along the axes, the B2 g orientation probes Cooper pairs where the nodes 
of the gap are located and does not sample regions where the gap is maximal. The 
reverse is true for the B\ g orientation. Thus the B2 g spectrum will reflect a smaller 
energy needed by the incoming light to break a Cooper pair and B\ g will have 
a larger energy scale. Lastly, since A\ g involves more of an average around the 
Brillouin Zone the energy scale can be less than that for B\ g . Moreover, screening 
effects for the fully symmetric charge distributions (^4i g ) can drastically reorganize 
the spectrum for this channel as the long range Coulomb interaction is brought into 
play. 

In most cases, the experimental Raman spectra for optimally hole-doped cuprates 
can be well modeled using an energy gap of d x 2_ y 2 symmetry. The relative peak po- 
sitions, low frequency power-laws and temperature dependence of the nearly elastic 
contribution for B\ g and B2 g scattering channels can be naturally explained assum- 
ing a d x 2_ y 2 paired superconducting state. The telling signature here is the presence 
of an uj 3 frequency dependence of the B lg channel frequencies compared to a linear 
dependence in other channels. This uniquely identifies the nodal positions of the gap 
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to lie along 45° in the C11O2 plane. However, it is to be noted that a small linear in 
frequency contribution is seen in the B\ g channcl in opt im al ly doped YBa2Cu307_^ 
(Y-123) which may be due to orthorhombic distortionsO which brcak the B\ g — A\ g 
symmetry distinction and/or may be due to the large Fano distortion of the back- 
ground due to the 340 cm -1 out-of-phase oxygen vibration£3 For the A\ g channel 
(which can not be purely selected in one particular geometry) the strong peak of the 
observed spectrum can still be fit by an appropriate choice of the gap and scattering 
amplitude but the theoretical prediction was found to be extremely sensitive to the 
numbcr of harmonics used to represent the k-dependence and thus to small changes 
in bând structure and/or small dopings. Moreover screening necessarily leads in 
the calculations to a much smaller signal than seen in experiments. It may be that 
additional physics is needed to describe the A\ g peak shape. However, the linear 
frequency behavior at small frequencies follows naturally due to a gap with nodcs. 
As in most corrclation functions for clean superconductors, only the absolute mag- 
nitude of the energy gap and not the phase is measured in Raman scattering. The 
current agreement of the calculations for clean d x 2_ y 2 superconductors compared 
to the data on optimally doped Bi2Sr2CaCu20s+5 (Bi-2212) are presented in Fig. 
1. 
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1. Comparison of the calculation: 
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performed for Raman scattering in clean d x 2. 



conductors compared to the data taken on optimally doped Bi-2212 by Staufer et al. 
5. 



y 2 super- 
in Ref. 



However, the theory cannot link a description valid for all temperatures to the 
normal state. For vanishing momentum transfers of the light, phase space for elec- 
tronic Raman scattering can only be opened due to scattering via impurities or 
electron-electron or electron-phonon interactions. In the superconducting state this 
phase space restriction is lifted since breaking Cooper pairs requires no net momen- 
tum transfer. 

Building upon the well known results of impurity effects in anisotropic s-wave 
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superconductors,tj Borkowski and flirschfeldla have also shown that while gap 
anisotropy is unaffected by impurity scattering in d— wave superconductors, sys- 
tematic impurity doping in anisotropic s— wave superconductors leads to thermally 
activated behavior as the gap anisotropy is averaged out and can thus provide an 
indirect way of determining if the gap has accidental nodal points or if the nodes are 
enforced by symmctry. Since smearing effects can obscure thermally activated be- 
havior, an observation of a threshold in energy, such as the density of states (DOS), 
could also be uscful. However, it is well known that DOS mcasurements can be 
problematic due to surface conditions and small coherence lengths. 

Moreover, many open questions are encountered in cuprate superconductors for 
the Raman spectra away from optimal dopingcâ In overdoped materials a small 
region (if any) of cubic behavior is scen in the B\ g channel before crossing over 
to a linear low frequency dependence. Strong scattering (~ uj) is seen down to 
the lowest frequencies measured in all overdoped materials for all other channels. 
In underdoped materials, dimculties are encountered to observe features associated 
with superconductivity in the B\ g channel while strong scatteringjs still seen at 
low frequencies for all channels and is linear in frequency as well.li3 Additionally 
the loss of overall scattering intensity of the B\ g channel relative tp_pther channels 
has been taken as evidence for a pseudogap in the normal stateJ13 This may be 
explained as the effect of losing part of the Fermi surface near (ir, 0) as hole pockets 
near (7r/2,7r/2) are developing. In addition, very distinct features and behaviors 
are seen in resonant Raman scattering when the incoming laser light is varied, and 
behavior reminiscent of two-magnon scattering in an antiferromagnetininsulator is 
seen to persist well into the underdoped region of superconductivity.liH Thcrcforc, 
while the existence of spectral weight down to the lowest frequencies implies the 
cxistcnce of a minimum energy gap at most as large as the energy resolution of 
the experiments, the strong symmetry dependence of the spectra as a function of 
doping can contain clues to the behavior of excitational dynamics and ultimately a 
mechanism of pairing for the superconducting state. 

While a complete theory now exists for Raman scattering in clean unconven- 
tional superconductors, this paper is devoted to augmenting the theory to include 
the effects of impurity scattering to exploit the fact that the response of a material 
to deliberate disorder can yield more information concerning the anisotropy of the 
energy gap but its phase as well. In this paper we will review the theory of Raman 
scattering in unconventional disordered superconductors with the goal of addressing 
how the scattering cross section in the superconducting state can provide informa- 
tion on the ground state symmetry as a function of doping. Open questions will be 
addressed and comparison to the available data will be made. 

This paper is outlined as follows: In Section 2 we shall review the essentials of 
the theory of Raman scattering in metals and confine our attention to the case of 
only one bând crossing the Fermi level. The relation of the symmetry of the Raman 
vertex and its coupling to an anisotropic energy gap is discussed. 

In Section 3 we will present a theory of Raman scattering in disordered uncon- 
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ventional superconductors. Isotropic s— wave impurity scattering is included in a 
gauge invariant way for superconductors with arbitrary gap symmetry at arbitrary 
temperatures. 

Section 4 presents the channel dependent spectra obtained for a d x 2_ y i and an 
anisotropic s— wave superconductor neglecting vertex corrections and are compared 
and contrasted. The symmetry aspects of the calculations as a function of fre- 
quency and temperature are discussed in detail. We also discuss the role of vertex 
corrections. In particular the existence of disorder generated collective modes are 
examined. 

Lastly, Section 6 compares the results of the calculations to spectra obtained on 
the cuprates at various dopings. The current agreement and lack of agreement of 
the data with the theory are discussed, and open questions are presented. A brief 
account of this work has appeared in Ref. 18. 

2. General formalism for metals 

Electronic Raman scattering measures effective charge fluctuations around the 
Fermi surface, 

P = ^7(k,q = 0) C j cCT c k:(T (1) 

k,CT 

with a scattering amplitude given by the Raman vertex, 

7 (k,q=0;w / ,w s ) = e I -e s + ! V (2) 

m * — ' 

(n, k | e s • p | v, k) {v, k | e 1 • p | n, k) (n, k | e 1 • p | v, k) (v, k | e s • p | n, k) " 



e(k) - e^(k) + ui e(k) - e I/ (k) - lj s 

where e , e s (uj], ujs) are the incident and scattered photon polarization vectors (en- 
ergies), p = —iKV, and e(k). and e„(k) are the Bloch conduction and intermediate 
state energies, respectivelyB While in the limit of vanishing light frequencies and 
for a single bând near the Fermi level the Raman scattering amplitude can be re- 
lated to the curvature of the conduction bând, in general this relation does not hold 
and is of questionable use for the cuprates. A thorough discussion of the Raman 
vertex is given in the Appendix of Ref. 10. An alternative approach is based on the 
experimental observation that the spectra near optimal doping in the normal state 
depends only mildly on the incoming laser light. Then the scattering amplitude 
can be taken as roughly independent of frequency and symmetry considerations 
can be given to its k— dependence. This of course misses resonant Raman scatter- 
ing which will be more relevant for smaller dopings nearer to the antiferromagnetic 
phase. However, it is a much more complex and unresolved problem of how to 
bridge non-resonant to resonant scattering. 

Therefore we elect to describe the Raman vertex 7 in terms of an expansion 
of Fermi surface (FS) or Brillouin zone (BZ) harmonics. The vertices depend on 
momentum throughout the BZ as: 

Big : 7(k) ~ cos(fca;a) — cos(k y a) + ■ ■ ■ (3) 
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i?2g ! 7(k) ~ sin(k x a) sm(k y a) + • • • 

A\ g : 7(k) ~ constant + cos(k x a) + cos(k y a) + • • • , 

for a 2-D lattice with lattice constant a. Here • • • represent higher order terms in the 
BZ expansion. This allows us to correctly classify the anisotropy and transformation 
properties of the scattering amplitude but leaves open the question of absolute 
intensities due to the unknown momentum independent prefactors of the expansion. 
The prefactors can be taken as parameters to fit overall intensities. Apart from A\ g 
(see Ref. 10) the prefactors have only a small effect on the frequency lineshape of 
the spectra. 

The channel-dependent Raman cross section is related to the channel-dependent 
Raman susceptibility x 7i7 via the fluctuation-dissipation theorem, 



d 2 <r LUs > 



dujdQ w\ 



r o S yn (q,u), 



S^ y (q,u) = [l + n(w)]Imx7, 7 (q,w), (4) 

7T 



with 



X 7 , 7 M - / ' dre-^(T T [p(r)p}), (5) 

Jo 

with T T the time-ordering operator and the imaginary part is obtained by analytic 

continuation, iu — > u + iO. Here ro = e 2 /mc 2 is the Thompson radius and we have 

set h = Ub = 1- 

In the absence of impurity scattering, the Raman response function can be 

written in terms of the k— dependent Tsuneto function 



w, x A ( k ) 2 , \E(k) 
A(k,î^) = -^^tanh 



2T 



1 1 



2E(k) + iu 2E(k) - iu 



(6) 



as 



x 7 ,iH = L r(k)A(k,w), (7) 



k 



with the excitation energy E 2 (k) — £, 2 (k) + A 2 (k), conduction bând ^(k) = e(k) — 
/i, /i the chemical potenţial, and energy gap A(k). We have neglected Coulomb 
screening (important for fully symmetric charge fluctuations A\ g ) as well as pair 
interactions responsible for maintaining gauge invariance and collective modea. 
details of calculations performed with Eq. (uh, the reader is referred to Refs. QH 

From Eqs. ([w]) , we see that the Raman vertex couples directly to the energy gap 
under the momentum summation. Since the B\ g channel assigns maximum weight 
along the BZ axes (0, ±1) and (±1, 0), while B^g weights the diagonals (±1, ±1), the 
orientation of the incident and scattering light polarizations thus choose effective 
charge fluctuations on the corresponding regions of the FS. Therefore, if for instance 
a d x 2_ y 2 energy gap is used, A(k) = A [cos(k x a) — cos(k y a)]/2, the B lg {B2 g ) chan- 
nel samples regions of the gap maximum (minimum) , respectively. The combination 
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of the two symmetries thus gives informaticul about the nodal behavior as well as 
the maximum value of the energy gap Ao- 

3. Disordered unconventional superconductors 

3.1. T-matrix approach 

We now consider the effect of scattering by spinless, noninteracting, isotropic 
impurities on the Raman susceptibility of unconventional superconductors. We use 
the seif consistent T— matrix approach to incorporate repeated scattering events 
from a single impurity site and dress both the Green's functions as well as the 
vertex. 

The two parameters in this theory are the cotangent of the scattering phase shift, 
c = cot (<5), and the impurity concentration n, described through the scattering rate 
r = (N/V)rii/irNF, where N/V is the electron density and Np is the density of 
states per spin at the Fermi level.t^l These enter into the T-matrix and the seif 
energy in particle-hole (Nambu) space: 

t(k,iu) = TT(k,k,iu), (8) 

The T matrix (in terms of the seif energy) is depicted in Fig. 2 and satisfies a 
Bethe-Salpeter equation, 

f (k, p, iu) = V(k, p)+Y, V(K k')G(k', iu)f (k', p, iu), (9) 

k' 

with G the Green's function in Nambu space which contains the T matrix via 
Dyson's equation. The matrix V(k, p) is the impurity scattering potenţial for a sin- 
gle electron, taken to be independent of the electron's spin. Moreover, the potenţial 
is taken to be independent of momentum as well to represent s—wave scattering 
only. Therefore the seif energy is k— independent as well. 

Expanding the seif energy in quarternions, S = y^-_ Q TjEj, where fo is the 
2x2 unit matrix and fi(i — 1, 2, 3) are the Pauli matrices, the one-particle Green's 
function can be written as 

G (k)H = !^±ISÎ3±ÂWn_ (10) 

(^) 2 -£ 2 (k)-A 2 (k) 

The tilde indicates the renormalized frequency, gap, and bând energy via 

iu = iu-'Eo(iu), Ă(k) = A(k) + S 1 (w), f(k) =£(k) - E 3 (w). (11) 

The matrix seif energy is given in terms of the integrated Green's function gi (iu) = 
2^E k Tr{nG(k,iu)} as 



, 9 o(iu)To - gi(iu)ri - ct 3 
c 2 - gl(iu) + gl(iu) 



Y>(iu) = T — 2 A — =EotoH-EitiH-E 3 t3. (12) 
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T = 



Fig. 2. Diagrams for the T-matrix seif energy £ and the renormalized vertex V as defined in the 

text. 



Here Tr denotes taking the trace. The off-diagonal seif energy Ei is zero only for 
odd-parity states or states which possess reflection symmetry such as d x 2_ y 2 or 
d xy . Here we have assumed particle-hole symmetry and further set gs(iu;) = 0. 
This allows us to cast the seif energy in the simplified form given in Eq. (n2h . As 
discussed in Ref. 19, this approximation is valid for either weak (Born, c >> 0) 
or strong (unitary, c = 0) scattering provided that impurity vertex corrections play 
only a minor role. This is shown to be the case in Section 4.3, and therefore we use 
Eq. (O) for the seif energy. For a more detailed discussion of this point, the reader 
is referred to Ref. 19. 

The expression for the Raman response in the limit of vanishing momentum 
transfer is given by 

X 7 , 7 (q -> 0, tîî) = -TJ2 E T>{ 7 (k)f 3 G(k, ^ + )^ 3 f (k, iCo+, iu>-)G(k, iu;-)}, 

iu k 

(13) 
where ui = ui ± fî/2, and 7T3,T are the bare and impurity renormalized Raman 
vertices, respectively. The diagrammatic representation for the renormalized vertex 
is shown in Fig. 2 and can be expressed as the integral equation 

f(k, iu> + ,iQ-) = 7(k)f 3 + Y^ T + G{p, iu> + )î(p, iui + ,iu}-)G(p, iu)-)f-, (14) 



with f ± = f (iui*). 

Eqs. (8-14) form a closed set of equations for the Raman response of unconven- 
tional superconductors in the seif consistent T-matrix approximation. We remark 
that while the impurities are included in a gauge invariant way, the neglect of the 
pairing interaction in the renormalized vertex does not satisfy gauge invariance 
and therefore all information regarding the existence of pairing interaction induced 
collective modes and sum rules is lost. While a gauge invariant treatment is pos- 
sible for disordered s-wave superconductorsu a similar treatment has not yet been 
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performcd for unconvcntional superconductors. We note that for clean supercon- 
ductors, the gauge invariant Raman response has been calculated in Ref. 9 for 
d~ wave superconductors. There it was found that the collective modes which do 
exist in different Raman channels are damped and lead only to small shifts in the 
relative peak positions of the Raman spectra in each channel. The low frequency 
behavior in particular is unaffected. However the Goldstone mode which appears as 
a consequence of the spontaneously broken U(l) gauge symmetry, provides for both 
longitudinal and transverse screening of the Ai g response in a non-trivial way, as 
discussed in length in Ref. 9. Therefore we are not in a position to discuss the be- 
havior of the A\ g Raman response in the presence of impurities for unconvcntional 
superconductors since this would require a gauge invariant treatment. 

3.2. Solution of the integral equations 

To solve Equation (|14|), it is convenient to first remove the k— dependent term 
and define T(k, iu), iuo') — f 3 7(k) + 5{iu. iiv'), and then expand once again in quar- 
ternions, 6 = ^2 i=0 TiSi. Eq. ( |14| ) thcn turns into a 4 x 4 matrix integral equation 
with no non-zero elements for general A(k). To simplify a solution, we now restrict 
attention to odd-parity states or states which possess reflection symmetry such as 
d x 2_ y 2 or d xy . Therefore the energy gap satisfies Ş^ k A(k) = and the off-diagonal 
term of the seif energy Ei can be set to zero. The resulting matrix equation then 
simplifies to 

S = J2lIooL 00 (k)5 + 7 03 (L 33 (k) 7 (k) + L 33 (k)S 3 )} 

k 

*i = ^[hiL 11 ^! - 7 12 (L 23 (k) 7 (k) + L 22 (k)S 2 )] 

k 

<5 2 = ^[/ii(i 23 (k) 7 (k)+i 22 (k)^) + / 12 L 11 (k) ( 5 1 ] 

k 

S 3 = 5>oo£ 33 (k)( 7 (k) + 6 3 ) + I 03 L 00 (k)5 }, (15) 

k 

where we have simplified notation and defined the following: 

Ţ rp + rp— | rp+rp— j T^+T^ - rp-\-rp — 

Lm - J-n J-Q + 1 3 1 3 , Jn — i i - 1 3 1 3 , 

% = T+T 3 - + T+T - , I 12 = -i{T+T^ - T+T Q - ) , (16) 



and 



L 00 (k) = 



llr v IUI' IUI 



+ A 2 (k)+e+(k)r(k) 
iV+(k)JV-(k) 

+ A 2 (k)-|+(k)r(k) 



L 22 (k) 



N+(k)N-(k) 
iCj+iCu- - A 2 (k) - |+(k)r(k) 



N+(k)N-(k) 
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33 _ ^ + ^--A 2 (k)+|+(k)|~(k) 

W JV+(k)iV-(k) 

23 _ iAQL)[iu+ - iu,'] 
L W ~ AT+(k)iV-(k) ' (1?) 

with AT±(k) = (tw*) 3 - (^(k)) 2 - A 2 (k). 

The matrix equation can then be solved for general vertex 7(k). Dcnoting 
(j4(k)) = X^k^C 4 )' we obtain the solution to trie vertex equation as 

< 7 (k)L 33 (k)) +( 7 3 (L 33 (k)} 

00 — -<03 — 



Si 



1 - /oo(£°°(k)) 

( 7 (k)£ 23 (k))+£ 2 (L 22 (k)) 
l-J u (L"(k)) 



5 2 = < 7 (k)L 23 (k)> " i-'"<i"W) 



1 r,/r22n : ^ i 7-2 (L"(k)>(L"(k)) ' 
l-iii^L (kj; +7 12 îzT^irrp^ 

r /n 2 ,{£°°(k)) 

U l /oo(L 33 (k)) - /gs ^gffg 

The conserving approximation for the vertex equation (|14| ) cnforces that the impu- 
rities have been treated in a gauge invariant way for the the renormalized vertex. 
However, the full theory is not gauge invariant since the pairing interactions respon- 
sible for superconductivity have not been included in the vertex equation. 

3.3. Raman susceptibility 

We are now in a position to obtain our final result for the Raman susceptibil- 
ity. Substituting Eq. ( |l8| ) into Eq. ( |13| ) and performing the trace, the Raman 
susceptibility follows as 



^q = 0,zfi) = 2T]TJ( 7 2 (k)L 33 (k)) 



-< 7 (k)L 33 (k)> 



2 


T | r2 <i°°(k)) 


1- 

2 


j /rmnA\ n <£ 00 (k)}<L33 (k )) 

j | r2 (^"(k)) 
J ll + i 12l-/ 11 (LH(k)) 



-( 7( k)L 23 (k)) 2 ^ - ,;_— ^y;;;^,^ >. ai» 

L (k)) 



l-In(^(k))+I? 2 ™S> 



The first term is the response in the absence of vertex corrections while the next 
two terms are due to vertex corrections. To obtain the final result for the Raman 
cross section, one must analytically continue by letting iQ — > Î7 + iO and take the 
imaginary part of Eq. ( |19|) to be put into Eq. <Q). 

We now consider the role of symmetry in the vertex corrected response. We focus 
on Ai g scattering first to show its connection to gauge invariance. For isotropic 
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density fluctuations, 7(k) = constant, and thus by symmetry the last terni in Eq. 
( |l9| ) is zero. Part of the second terni is canceled by the first term and we are still 
left with a finite response at q = for isotropic density fluctuations. This is due 
to having a theory which breaks gauge invariance and violates the /—suni rule and 
thus partide number conservation. The density response will be made to vanish 
when the pair interaction is included in the vertex renormalization and when long- 
range screening by the Coulomb forces are taken into account. This only affects the 
results for fully symmetric scattering (Ai g ). We note that for general scattering in 
the A\ g channel, 7(k) need not be a constant. However, by symmetry we see that 
the third term in Eq. (|19|) vanishes for A\ g scattering in d— wave superconductors. 

Next we consider B\ g and B 2g scattering channels. Again by symmetry we 
see that the second term vanishes for both channels. However, the third term 
contributes for channels with the same symmetry as the energy gap A(k). For a 
gap with d x 2_ y 2 pairing, the third term is zero for B 2g scattering but contributes 
for Big. The opposite is true for d xy pairing. Thereby we see that the vertex 
corrections are very symmetry dependent and can by themselves lead to channel- 
dependent line-shapes. 

In the following section we will discuss the results of the theory in the absence of 
vertex corrections and at zero temperature to bring out the differences between the 
results from clean unconventional superconductors. In particular we will investigate 
how impurity effects can help determine the phase of the energy gap. We defer 
discussions of the role of vertex corrections to the following section. 

4. Scattering in the absence of vertex corrections 

4.1. T=0 results 

In this section we will calculate the effects of impurities on the Raman spectra 
for unconventional superconductors by evaluating Eqs. (|l7r~|l9|). To simplify the 
calculations, in what follows we will work with an approximately cylindrical FS and 
restrict all momentum averages to the FS. The vertices then depend on azimuthal 
angle <\> around the FS as 

B lg : 7 (0) = V2cos(20) (20) 

B 2g : 7 (0) = \/2sm(20), 

where we also neglect higher order terms. We also neglect bând structure details 
and for simplicity take an infinite bând as well. We thus approximate all k sums as 
an energy integral times an angular integral around the FS: 



Ş^-> n f I dţdn^. 



The reader is referred to Ref. 9 for a discussion of consequences of the these ap- 
proximations. For instance we are no longer able to discuss the role of the van Hove 
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singularity. Recent calculations show that the van Hove singularity due to the flat 
bands located near (tt, 0) shows up as a peak in the B\ g Raman spectrum at an 
energy ~ -y/A 2 + e^ H , where e v u is the distance in energy of the flat bând from the 
Fermi level near the zone axis (see Ref. 10 and Branch and Carbotte in Ref. 8). 
However no such van Hove peak, which should be sensitive to bând structure and 
should not be tied to the energy gap, is seen in experiments. Since the van Hove is 
at higher energies for most systems it is likely that quasiparticle inelastic scattering 
is sufflciently strong to damp this feature if it exists, or it may be lying at furthcr 
distances away from the Fermi level and thus have inappreciable residue. In any 
case we neglect this since impurities are sufficient to suppress the van Hove feature 
in any case. 

In this section we consider specifically a d x 2_ y i paired superconductor, with 
energy gap Ad(</>) — Ao cos(20) and compare with the results for a hypothetical 
anisotropic s— wave superconductor with A s ((/>) = Ao | cos(20) |. In the absence 
of impurity scattering, the channel-dependent Raman spectra (as well as the re- 
sults for other correlation functions) would be indistinguishable. We will consider 
only the T = response and defer consideration of finite temperatures as well as 
vertex corrections to the next subsection and following section, respectively. Since 
Ş^ k A(k) is not zero for the anisotropic s case, we reinsert the Si term into L 33 (A 
is replaced by A in Eq. (|17|)) when evaluating integrals for the vertex uncorrected 
calculations. This is the only modification needed. 

The results of the calculations for unitary scattering c = are shown in Figs. 
3 and 4 for the B\ g and Big channels for both superconductors. One immediately 
notices that the impurities smear out the logarithmic divergence of the peak at 2Ao 
for the B\g channel for both superconductors. Moreover, as the scattering rate 
r/Ao increases past roughly 0.5, the diffcrcnccs between the peak of the spectra 
in both channels become less pronounced as the scattering becomes the dominant 
energy scale in the problem. But as a consequence of the gap renormalization for 
anisotropic s-wave superconductors, Si is non zero and the gap becomes averaged 
out by the disorder and a threshold develops at uj g = 2A m i„. For small impurity 
scattering, this leads to a reduction of the relative peak positions for the B\ s and 
B2 g channels compared to d x 2_ y 2, as shown in Figs. 3 and 4. As the disorder 
is increased, the peak positions will coalesce and the spectra recover a channel 
independent, isotropic s— wave form, □ 

X!'-™««M ~ e (^ - 2A aue )^p for u ~ 2A avc 

with A ave = 2A /7T. 

However, the main difference for d x i_ y i and anisotropic s-pairing lies in the low 
frequency behavior, which is very channel dependent, as shown in the insets of Figs. 
3 and 4. For clean materials, both superconductors show the characteristic w(o; 3 ) 
behavior for the B 2 g(Bi g ) channels, respectively. 

For d x 2_ y 2 pairing, while impurities do not change the linear in frequency behav- 
ior for the B^ channel, below a characteristic frequency lo* ~ VTA the behavior 
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Fig. 3. Raman spectra for B\ g and i?2g channels for d x 2_ y 2 pairing for T = 0. Values of T/A 
are shown for resonant scattering. The inset shows the low frequency behavior on a log-log scale. 
The Big spectra show a cross over from ui to ui 3 at roughly ui* ~ vTA. The bounding lincs in 
thc Big insets are guidcs for ui (upper) and a; 3 (lowcr) behavior. 
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Fig. 4. Raman spectra for B\ g and i?2g channels for a supcrconductor with a | d x 2_ y 2 | gap 
symmetry. Inset shows the low frcquency behavior which is dominatcd by a threshold for both 
channels. 
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ii g channcl. This is due to a nonzero density of 
statesLat the Fermi level, which allows for normal-state-like behavior to be recov- 
ered. EJ As in the case of the penetration depth,EJ the scale w* grows with increasing 
impurity concentrations, and will be shown in the next section to be strongly tem- 
perature dependent. However, the exponent is symmetry dependent (remains 1 for 
B2 g and A\ g channels, while decreases from 3 to 1 for the B\ g channel). This is 
in marked contrast to the impurity dependence of the spectra for anisotropic s gap 
as seen in Fig. 4. The low frequency behavior is dominated by the threshold in 
this case for all channels. Moreover, the impurity dependence of the B\ g chan- 
nel is opposite to what one would expect if the gap was anisotropic s-wave. For 
anisotropic s, one sees a transfer of spectral weight out to higher frequencies due to 
the development of the threshold while for d x i_ y i the transfer of spectral weight is 
towards lower frequencies. Thus impurity scattering can provide a clear qualitative 
way of distinguishing energy gaps of different symmetry. 

It can be shown that Born scattering (c >> 0) leads to the same effects except 
for low frequencies. The spectra reorganize at higher frequencies (u ~ 2Ao) in the 
same way as for resonant scattering. Also, as has been pointed outEJ there is little 
difference between Born and unitary scattering for anisotropic s— pairing even at 
low frequencies. However, the low frequency behavior is changed for d— pairing. In 
the Born limit, the low frequency behavior in the Bi g channel crosses over from u) 
to w 3 at a much smaller frequency to* . Therefore the low frequency behavior of the 
spectra is similar to the behavior obtained in the absence of impurities except for 
only very small frequencies in the B\ g channel. A substanţial linear behavior for 
Bi g can be obtained in the Born limit for large scattering, T/Ao ~ 1 + c 2 , but such 
a large scattering almost completely smears out the peak feature at 2Ao- 

For small scattering, the threshold generated in anisotropic s— wave supercon- 
ductors can be obscured via inelastic scattering or experimental resolution. How- 
ever, the "effective" exponents for all channels would grow as the threshold develops 
for increased impurity scattering. Since low temperatures and low frequencies can 
be achieved at high resolution (~ 8 cm -1 ), the smearing can be controlled to allow 
for a systematic check of the impurity dependence to determine whether the gap 
has accidental or intrinsic zeroes. 

4.2. Finite temperatures 

We now consider the effect of finite temperatures on the calculations. The same 
physics that is manifest in the symmetry dependent low frequency behavior gives 
rise to a similar characteristic temperature dependence which can also be used to 
ascertain the energy gap symmetry. 

We first note that for temperatures above T c , the Raman spectrum for channcl 
L is given by a simple Lorentzian: 

o T 
X£(n, T > T c ) = 2iV F7 | ( ^ )2 + i , (21) 
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where r = 1/2T is the impurity scattering lifetimeE3 In the presence of vertex 
corrections r becomes channel dependent. This spectrum peaks at frequencies Î7 = 
I/t, rises linearly with frequency and falls off at large frequencies as 1/0. More 
complicated forms for r can be obtained for k— dependent scattering and when 
electron-electron interactions are taken into account.tJ 

Our results recover this form as T approaches T c . Therefore, any cubic rise of 
the B\g channel must shrink as temperatures increase and the peak in the spectra 
must approach the normal state value as well. Both factors lead to a temperature 
dependence of the crossover frequency uj* estimated in the previous subsection. 



2N C 




O 0.5 



co/A 



Fig. 5. -Bi 9 Raman spectrum evaluated for T/A = 0.05 and for unitary scatterers at successive 
temperatures T/T c = 0.1,0.5,0.8,0.9,0.99 and 1, respectively. Here a weak coupling temperature 
dependence of the energy gap A/T c = 2.14c! has been assumed for simplicity. 



The temperature dependence of the results for d x 2_ y 2 superconductors are shown 
in Figs. 5 and 6. Here we have assumed a weak coupling temperature dependent 
energy gap for simplicity, whose weak coupling value is 2Aq/T c — 4.28. From the 
Figurc it can be seen how the normal state Lorentzian line-shape rises out of the su- 
perconducting response as the temperature increases and the energy gap decreases. 
An additional peak appears at ii = 1/t in both channels with little weight at low 
temperatures and increasing in magnitude as the temperature increases. The low 
frequency behavior is thus dramatically changed so that the spectrum can easily 
loose any trace of cubic behavior if the impurity scattering and/or the temperature 
is large enough. Therefore only at very low temperatures can an elastic scattering 
rate be inferred from the data. 

Additional information can be obtained by observing the low frequency behavior 
as a function of temperature. The ratio of the superconducting response to the 
normal state response in the static limit SI — v shows characteristic temperature 
power-law behavior in the case of clean d— wave superconductors as discussed in 
Refs. 7 and 9. There it was shown that for clean materials the ratio is given by a 
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Fig. 6. B2 g temperature dependent spectrum using the same parameters as in Fig. 5. 



weighted average over the FS of a Fermi function / as 

*?...(« ->0,T) 2(/(| A(k) |) | 7 (k) | 2 ) 



x'ijn -> o,t) 



(I 7(k) | 2 ) 
T 3 , T^O, 
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B 



iaj 



T^O, 



B 



2y 



(22) 



where the last two lines gives the same exponents for the low temperature behavior 
as the exponents for the frequency dependence. It can be seen that the spectra 
develop a suppression of states at low frequencies faster for the B\ g channel than 
B2 g , which has been borne out in experiments on clean cuprate samples.El 

In the presence of disorder this simple expression does not hold. This is due 
to the fact that the Green's functions do not have an undamped simple pole any 
longer due to the pair breaking nature of impurities in a d— wave superconductor. 
Moreover, for clean systems above T c there is no Raman scattering so this limit is in 
a sense artificial. The inclusion of impurity scattering verifies the above form if the 
limit of vanishing scattering is taken as to make the definition sensible. However, 
for finite scattering the static limit ratio is not given by Eq. (|2l|). For d x 2_ y i 
pairing, since for all temperatures below T c the Big and A\ g channels have a linear 
dependence on fi for small frequencies while even B\ g has a small linear dependence 
as well, this linear dependence mimics the normal state behavior and thus gives a 
constant term even at T = 0. Therefore all curves will be shifted upwards at low 
temperatures. 

The static ratio for various values of disorder is shown in Figs. 7 and 8 for 
unitary scattering for a d x 2_ y i paired superconductor. The results are compared to 
data taken on optimally doped Bi-2212. The number this ratio takes for T = can 
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Fig. 7. Big static ratio as defined in the text as a function of progressively stronger disorder 
(r/A(T = 0) = 0,0.01,0.05,0.1,0.5 and 1 from bottom to top). The circles are takcn from the 
data of Hackl et al. in Ref. 15. 



be calculated analytically. If we define 6 as 5 = £q (w — ► 0) the ratio becomes 

7 2 (k) 



x?. a .(n->o,r = o) . 3 ( ^ţaw^ ) 



x£. a .(n->o,r = o) 



<7 2 (k)) 



(23) 



This relation shows that in the limit of vanishing scattering, S << Ao, once again 
we recover channel dependent exponents. For the B\ g channel, Eq. (J23|) vanishcs 
as (i5/Ao) 3 while for B^g it vanishes linear ly in 6. 

Finally we determine how large S is for both Born and unitary scattering by 
solving Dyson's equation seif consistently.O For the case of Born scattering c >> 1, 
then 6 ~ e °' r and is exponentially small. For unitary scattering c = 0, S is 
given as a seif consistent solution of S ~ e 5 ", which is the Omega function W 
S = 2 v /r/A /\/— 2W(— r/A ). The VF function provides logarithmic factors and 
therefore 5 ~ yT/Ao in the unitary limit. Therefore for unitary scattering the 
constant contribution to the static response can be quite large and therefore gives 
a Raman spectrum which produces normal state behavior with an albeit reduced 
intensity even at T = 0. We note additionally that a constant contribution to the 
ratio of 20% and 33% has been seen for the B\ g and Bi g channels, respectively, in 
Bi-2212 by Yamanaka in Ref. 5. 

In closing this section we remark that all calculations have been performed 
using a weak coupling temperature dependence of the energy gap, which does not 
provide a good description for the cuprates. However it is trivial to generalize these 
results to strong coupling energy gaps by using the interpolation formula provided 
in Appendix B of Ref. 9, which gives the temperature dependence of the gap in 
analytic form using the specific heat jump and angular averages of the energy gap as 
input parameters. The above results are only quantitatively changed near T = T c . 
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Fig. 8. Big static ratio as defined in the text as a function of progressively stronger disorder 
(r/A(T = 0) = 0, 0.01, 0.05, 0.1, 0.5 and 1 from bottom to top). The circle is taken from the data 
of Hackl in Ref. 15. 



Moreover, we have performed all calculations for a 2-D cylinder-like FS and 
therefore do not capture log corrections that would occur for all low frequency and 
temperature behavior of all quantities in 3-D. In the following subsection we will 
discuss the role of vertex corrections which have been neglectcd in this subsection 
as well. 



4.3. Vertex corrections 

We return to the expression Eq. (09) for the vertex-corrected Raman response 
and discuss what changes to the spectra occur when impurity vertex corrections 
are taken into account. In particular we discuss whether any disorder generated 
collective modes appear in the spectra and examine any channel dependence which 
results from the vertex corrections. We remind the reader that we are neglecting 
vertex corrections from the particle-particle pairing interaction and therefore these 
results are not truly speaking gauge invariant even though the impurities themselves 
are treated gauge invariantly. 

We limit discussion only to the case for T = 0, where the Matsubara sum can 
be converted into an integral and the analytic continuation to the real frequency 
axis can be done without having to resort to Pade approximants or other numerical 
techniques such as maximum entropy. 

We first examine whether there exist any collective modes introduced by the 
disorder for a d x i_ y i superconductor. We note again that the second term in Eq. 
( |l9|) is non-zero only for A\ g scattering and therefore we only will focus on the third 
term which is non-zero only in the B\ g channel. Therefore there are no impurity 
vertex corrections for the B 2g channel and subsequently no collective modes. 

We have searched for frequencies where the real part of the denominator is zero 
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in the second term of Eq. (M) and have not been able to identify the existence of 
any collcctivc modcs. We plot in Fig. 9 the results of the vertex corrected response 
for the Big channel and compare them to the results we obtain neglecting vertex 
corrections. We see that the spectra are essentially unmodified apart from a sup- 
pression of spectral weight for frequencies near the gap edge 2Arj. In particular we 
do not see a well defined collective mode and the low and high frequency behaviors 
remain essentially unchanged. Therefore the neglect of vertex corrections does not 
lead to any substanţial changes to the impurity averaged Raman response. This 
is similar to the results obtaincd for the pair-interaction-vertex-corrected response 
obtained in Ref. 9, whcre it was shown that the vertex corrections lead to essentially 
no major modifications. 
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Fig. 9. Comparison of the calculation of the B\ g zero temperature spectrum with and without 
vertex corrections for T/A = 0.1 and unitary scatterers. 



5. Comparison to data on cuprate superconductors. 



5.1. Optimally and overdoped cuprates. 

In this section we compile our results and compare them to recent data taken 
on the high T c cuprates. We remark at the onset that no systematic check of the 
effects of dopant impurities on the electronic Ramam spectrajmve been undertaken 
in the metallic state. We do note that doping of Zn lill and PrEB have been inspected 
via Raman scattering on the two-magnon spectra obtained in insulating and under- 
doped Y-123. There it was seen that the two-magnon feature decreased with the 
introduction of Zn impurities, while a spectral reorganization of the Ai g intensity 
was observed at roughly a temperature 20 percent higher than T c . Since Zn is 
believed to seriously distort the local antiferromagnetic order in the CuC>2 plane, 
this finding lends support that antiferromagnetic correlations are strong in these 
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compounds for doping levels where superconductivity is establishcd. The cause of 
the spectral reorganization for Pr doping is unknown at present. 

Thcre is recent evidence which suggests that these materials are intrinsically 
disordcrcd in the overdoped regions of their phase diagram. Measurements of Ra- 
man scattering in the normal state of overdoped Bi-2212 by Hackl et al. in Ref. 
15 have infcrrcd a larger extrapolated zero temperature scattering rate than that 
obtained for optimal doped samples. Moreover the scatteri ng rate obtained from 
the Big channel matched that obtained from D.C. transporţi^ Also, recent muon 
spin rotationEU and Hali effect dataEi have suggested that at least for overdoped ma- 
terials there seems to be intrinsic disorder which is manifest in larger quasiparticle 
scattering rates. 

One possible scenario of how overdoped materials can be considered to be more 
disordered lies in a crossover argument from 2-D to 3-D behavior. In optimally 
doped materials, the carriers responsible for superconductivity are confined to the 
2-D CuC>2 layers. It is known that these material^become more isotropic (c-axis lat- 
tice parameter decreases) as they are overdoped.Ea This is manifestin that the c-axis 
conductivity becomes more and more metallic with overdoping£3 As the material 
becomes more 3-D the marginally confined carriers can interact with structural dis- 
tortions in the charge reservoirs which separate the CuC>2 layers. This has been put 
forth in Ref. 27 as an explanation as to why Ca doping Y-123 (which substitutes 
Y) leads to a much slower decrease in T c then Zn doping (which substitutes Cu). 

Thcrefore we will apply the theory as a possible explanation for the channel 
dependent Raman spectra in overdoped materials. We remark once again that we 
will treat the energy gap magnitude as a phenomenological parameter to be used to 
fit the poşition of the B\ g Raman peak and assume d x i_ y i pairing to be independent 
of dopingEil. This is not entirely unreasonable since a common feature of the Raman 
data at any polarization is the observance of large intensity at frequencies extending 
towards vanishing energy transfers, implying the existence of gap nodes. If one 
were to assume an isotropic s— wave energy gap, an explanation would be required 
in order to produce some smearing mechanism which is as large as 30% of the 
energy gap to produce the large scattering at low frequencies. Since the resolution 
of the Raman measurements is less than 10 cm -1 and inelastic scattering would be 
diminished by a factor of e~ A °/ T , this scenario seems unlikely. 

We present a fit of theory in Figs. 10 and 11 to the B\ g and B 2g spectra 
obtained on nearly optimally doped (T c = 86K), slightly overdoped (T c = 79 K), 
and appreciably overdoped (T c = 55 K) samples of Bi-2212 taken by Staufer et al. in 
Ref. 5 and by Hackl et al. in Ref. 15. No additional smearing mechanism is invoked. 
We see that while the peak positions of the data seem to move together as doping 
is increased the low frequency behavior of the combined data is well accounted for 
assuming a d— wave paired state and an increasing value of the resonant impurity 
scattering. This is manifest by the growth of the linear contribution of the B\ g - 
low-frequency response and the persistence of the B 2g linear fi behavior. Values 
of T/Aq estimated by the crossover frequency which separates linear from cubic 
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Fig. 10. Fit to the data taken for B\ g and i?2 9 channels in nearly optimally dopcd, as-grown Bi 
2212 (T c = 86 K, top panel) and slightly overdoped O2 annealed Bi 2212 (T c = 79 K, bottom 
panel) obtained by Staufer et al. in Ref. 5. Hcre r/A = 0.125(0.2) and A = 287(195)cm" 1 
have been used for the top (bottom) panel, respectively. Inset: Log-log plot of the low frequency 
portion of the B\ g response. The crossover frequency u;*/Ao = 0.38(0.45) for the top (bottom) 
panels, respectively. 
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bchavior in the B\ g channel yield uj* /A = 0.38(0.45, 0.58) for the optimally (mildly, 
appreciably over-) doped samplcs, rcspectively. Thus spectral wcight is bcing shiftcd 
to lower energies as the material becomes more overdoped. This lends supports to 
the conjecture that thcsc matcrials are more intrinsically disordcred for lower values 
of T c on the overdoped side of the phase diagram. However we note that the fit to 
the appreciably overdoped spectra in Fig. 11, especially the B 2g spectrum, is not 
as good as the fits to the other spectra. The position of the B 2g peak cannot be 
accounted for with our simple theory. Here additional physics seems to be needed. 
It may be possible that a different number of additional harmonics for the gap 
function (but still with d x 2_ y 2 symmetry) may move the peak outward (see Ref. 10 
and Branch and Carbotte in Ref. 8). However, without a detailed theory for the 
pairing mcchanism and its doping dcpcndcnce this remains an opcn question. 
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Fig. 11. Fit to the data taken for B\ g and i?2g channels in appreciably doped Bi 2212 (T c = 55 
K) obtaincd by R. Hackl et al. in Ref. 15. Inset: Log-log plot of the low frequency portion of the 
Big response. Here T/Ao = 0.22 and Aq = 120cm — 1 have been used. The crossover frequency 
ui*/A = 0.58. 



We have neglected quasiparticle inelastic scattering via e.g. spin fiuctuations or 
phonons which will lead to a furthcr smcaring of the curvcs and especially the large 
peak at 2 Ao in the B\ g channel. However, due to the rapid drop off of the scattering 
rates at low temperatures and frcquencies, the low frequency behavior of the Raman 
spectra will not be altered. Inelastic scattering, included phenomenologically in the 
calculations of Ref. 7 as well as in the calculations of Jiang and Carbotte in Ref. 8, 
is needed in order to reproduce the normal state behavior for temperatures above 
T c as well as the hat background seen at much higher cnergy shifts even at low 
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temperatures. Thus smaller values of disorder may then be used to fit the data. 
Moreover, the effect of resonant Raman scattering will be of more importance at 
higher frequencies and may play a pivotai role for a truly microscopic picture of 
Raman scattering as a function of doping. 

5.2. Conclusions and open problems. 

We close this section by summarizing our results and listing some open questions 
concerning a more complete picture of Raman scattering in the cuprate materials. 
We have seen that the theory can provide a satisfactory fit to the low frequency 
part of the channel dependent Raman response in the cuprate superconductors for 
optimal doping and for overdoping if we naively assume that the main effect of 
the doping is to introduce intrinsic resonant impurity scattering. This is certainly 
an open question and an unified picture of what happens even when materials are 
deliberately disordered via e.g. Zn or Ni doping is still lacking.tHl Moreover, the 
values of scattering needed to fit the Raman data are indeed quite large to describe 
the variation of T c as a function of doping if one applies Abrikosov and Gorkov's 
theoryS As noted, incorporated inelastic scattering could lead to smaller values of 
disorder needed to fit the data. However, the scattering rates implied by the Raman 
data match those obtained from D.C. transport, muon spin rotation, and Hali data 
rather well and do imply increased scattering on the overdoped side of the phase 
diagram. This still leaves the question of why T c is not zero unanswered. 

One possible scenario is that the scattering due to impurities is highly anisotropic. 
This evidence comes from the data as well, where both ţheRaman (see Stadlober et 
al. in Ref. 5 and Hackl et al. in Ref. 15) and Hali dataElllj suggest that scattering 
is larger along the BZ faces than along the diagonals. If this were true, part of the 
larger scattering could have the same symmetry as the energy gap and thus not 
be pair-breaking in the Anderson senseo The theory would have to be changed to 
incorporate an extended k— dependent impurity potentialEj 

This certainly appears to the case on the underdoped side of the phase diagram 
which we have not addressed. However, the scattering is most likely inelastic rather 



gring near "hot spots" 
ith underdoping as 
It is then clear that 



than elastic. Issues of anisotropic inelastic quasiparticle se 
of the FS have been invoked to explain the FS evolutio: 
well as transport rates obtained from a Boltzmann approach 
including inelastic quasiparticle scattering is a crucial missing point of the theory 
and is needed to explain the Raman spectra on the underdoped side of the phase 
diagram. In most Raman experiments, it is difficult to observe any superconductiv= 
ity related effects in the B\ g channel while they still persists for the B^ g responseJlj 
Moreover the relative intensities change as a function of doping in the normal state. 
For optimally and overdoped material, the Raman intensity in the B\ g channel is 
always larger than the Bi g intensity for any compound measured. However, the 
Bi g intensity has been observed to drop with underdoping and can become smaller 
than the B2 g intensity.E3 This may be related to the loss of the FS around the "hot 
spots" which are most effectively probed in B\ g polarization orientations via the 
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sclcction rules. Since Bi g measures the zone diagonals, the vestige of the FS, or 
"hole pockets" could still provide for a large electronic Raman signal. It is thus clear 
that the role of clcctron-electron interactions and incipient antifcrromagnetism will 
be needed to be incorporated into a theory of Raman scattering for underdoped 
materials. 

Moreover, recent measurementsEZI have shown that the Raman spectra for all 
channels in underdoped materials has a resonance pronie that is the same in the 
normal state as in the superconductor. Therefore one must include resonant Raman 
scattering processes into the theory even for underdoped materials. This is even 
more important if one wants to construct a theory for Raman scattering which can 
be extended to the insulating state and the two-magnon contribution. 

Thus in summary, a theory for Raman scattering in unconventional supercon- 
ductors can be useful to provide insight into the quasiparticle dynamics in high 
temperature superconductors for various regions of their phase diagram. However 
many issues are left unresolved and will require further work to build k— space 
anisotropies and electron correlations into the theory as well as resonant scattering 
processes. Application of the theory supports recent suggestions that these mate- 
rials are intrinsically disordered on the overdoped side of the phase diagram, and 
moreover d x 2_ y 2 pairing can provide an adequate description of the data. 

Acknowledgments 

T.P.D. would like to thank R. Hackl, G. Krug, M. Opel, R. Nemetschek, B. Stad- 
lober, J. C. Irwin, K. Hewitt, and J. Naeni for providing their data prior to pub- 
lication and for many useful and insightful discussions. T.P.D. would also like to 
thank J. C. Irwin and his colleagues at Simon Fraser University, where part of this 
work was completed. We would also like to acknowledge helpful discussions with 
D. Pines and A. Zawadowski. A.P.K. acknowledges support through a Heisenberg 
fcllowship of the Deutsche Forschungemeinschaft. 

References 



1. A. A. Abrikosov and L. A. Fal'kovskii, Zh. Eksp. Teor. Fiz. 40, 262 (1961) [Sov. Phys. 
JETP 13, 179 (1961)]. 

2. R. Sooryakumar and M. V. Klein, Phys. Rev. Lett. 45, 660 (1980). 

3. A. A. Abrikosov and V. M. Genkin, Zh. Eksp. Teor. Fiz. 65, 842 (1973) [Sov. Phys. 
JETP 38, 417 (1974)]; M. V. Klein and S. B. Dierker, Phys. Rev. B 29, 4976 (1984); 
H. Monien and A. Zawadowski, Phys. Rev. B 41, 8798 (1990). 

4. T. P. Devereaux, Phys. Rev. B 45, 12965 (1992); ibid. 47, 5230 (1993). 

5. S. L. Cooper et ai, Phys. Rev. B 37, 5920 (1988); R. Hackl et ai, ibid 38, 7133 
(1988); S. L. Cooper et ai, ibid., 11934 (1988); T. Staufer et al, Phys. Rev. Lett. 68, 
1069 (1992); A. Yamanaka et ai, Phys. Rev. B 46, 516 (1992); R. Nemetschek et ai, 
Phys. Rev. B 47, 3450 (1993); X. K. Chen et ai, Phys. Rev. Lett. 73, 3290 (1994); 
ibid. Journ. of Supercond. 8, 495 (1995); A. Hoffmann et ai, Physica C 235-240, 1897 
(1994); S. Donovan et ai, Journ. of Supercond. 8, 417 (1995); B. Stadlober et ai, J. 



9. 


T. 


0. 


T. 


.1. 


T. 


2. 


T. 


3. 


P. 


4. 


L. 


5. 


F. 



26 Raman Scattering in Cuprate Superconductors 

Phys. Chem. Solids 56, 1841 (1995); A. Yamanaka et al, in Spectroscopie Studies of 
Superconductors, eds. I. Bozovic and D. van der Marel, Proc. SPIE 2696, 276 (1996); 
A. Sacuto et al., LANL cond-mat preprint 9612052. 

6. J. Annett, N. Goldenfeld, and A. J. Leggett, to appear in Physical Properties of 
High Temperature Superconductors, Voi. 5, D. M. Ginsberg (ed.), (World Scientific, 
Singapore, 1996). 

7. T. P. Devereaux et. ai, Phys. Rev. Lett. 72, 396 (1994). 

8. X. K. Chon et ai, Physica C 227, 113 (1994); 227, 1089 (1995); M. C. Krantz and 
M. Cardona, J. Low Temp. Phys. 99, 205 (1995); D. Branch and J. P. Carbotte, Phys. 
Rev. B 52, 603 (1995); 54, 13288 (1996); C. Jiang and J. P. Carbotte, Solid State 
Commun. 95, 643 (1995); Phys. Rev. B 53, 11868 (1996); W. C. Wu and A. Griffin, 
Phys. Rev. B 51, 1190 (1995); 52, 7742 (1995); O. Narikiyo and K. Miyake, J. Phys. 
Soc. Jpn. 64, 3579 (1995); D. Einzel and R. Hackl, J. of Raman Spect. 27, 307 (1996); 
M. Cardona et. ai, in Spectroscopie Studies of Superconductors, eds. I. Bozovic and 
D. van der Marel, Proc. SPIE 2696, 182 (1996); T. P. Devereaux, ibid., p. 230. 

. P. Devereaux and D. Einzel, Phys. Rev. B 51, 16336 (1995); ibid. 54, 15547 (1996). 
. P. Devereaux, A. Virosztck and A. Zawadowski, Phys. Rev. B 54, 12523 (1996). 
. Strohm and M. Cardona, preprint. 

. P. Devereaux, A. Virosztck and A. Zawadowski, Phys. Rev. B 51, 505 (1995). 
Hohenberg, Zh. Eksp. Teor. Fiz. 45, 1208 (1963) [Sov. Phys. JETP 18, 834 (1964)]. 
S. Borkowski and P. J. Hirschfeld, Phys. Rev. B 49, 15404 (1994). 
Slakey et ai, Phys. Rev. B 42, 2643 (1990); X. K. Chcn et ai, Phys. Rev. B 48, 
10530 (1993); T. Katsufuji et al, Phys. Rev. B 48, 16131 (1993); R. Hackl et al, J. 
Low Temp. Phys. 105, 733 (1996); in Spectroscopie Studies of Superconductors, eds. 
I. Bozovic and D. van der Marel, Proc. SPIE 2696, 194 (1996); C. A. Kendziora, ibid., 
p. 223; Phys. Rev. B 52, 9867 (1995); Phys. Rev. Lett. 77, 727 (1996); K. Hewitt et 
al, to appear in Phys. Rev. Lett. 

16. J. C. Irwin et al. submitted to Phys. Rev. Lett.; R. Hackl et al. submitted to Phys. 
Rev. Lett. 

17. G. Blumberg et al, in Spectroscopie Studies of Superconductors, I. Bozovic and D. 
van der Marel (eds.), Proc. SPIE 2696, 205 (1996); M. Kang et al, Phys. Rev. Lett. 
77, 4434 (1996) and references therein. 

18. T. P. Devereaux, Phys. Rev. Lett. 74, 4313 (1995); Journ. of Supcrcond. 8, 421 (1995); 
J. Phys. Chem. Solids 56, 1711 (1995). 

19. P. J. Hirschfeld et al, Phys. Rev. B 37, 83 (1988). 

20. K. Ueda and T. M. Rice, in Theory of Heavy Fermions and Valence Fluctuations, 
edited by T. Kasuya and T. Saso (Springer, Berlin, 1985); L. P. Gorkov, Pis'ma Zh. 
Eksp. Teor. Fiz. 40, 351 (1994) [Sov. Phys. JETP 40, 1155 (1985)]. 

21. P. Hirschfeld and N. Goldenfeld, Phys. Rev. B 48, 4219 (1993); F. Gross et al, Z. 
Phys. B 64, 175 (1986). 

22. A. Zawadowski and M. Cardona, Phys. Rev. B 42, 8798 (1990). 

23. C. M. Varma et al, Phys. Rev. Lett. 63, 1068 (1989); B. S. Shastry and B. I. Shraiman, 
Int. J. Mod. Phys. B 5, 365 (1991); A. Virosztek and J. Ruvalds, Phys. Rev. Lett. 67, 
1657 (1991); V. N. Kostur, Z. Phys. B 89, 149 (1992); A. P. Kampf and W. Brenig, Z. 
Phys. B 89, 313 (1992). 

24. P. A. Lee, Phys. Rev. Lett. 71, 1887 (1993). 

25. M. Rubhausen et al, preprints. 

26. C. Kendziora et al, Phys. Rev. B 45, 13025 (1992); ibid. 48, 3531 (1993). 

27. C. Bernhard et al, Phys. Rev. Lett. 77, 2304 (1996). 

28. N. P. Ong in Physical Properties of High Temperature Superconductors, ed. by D. 
M. Ginsberg (World Scientific, Singapore, 1990), Voi. 2; J. P. Rice, J. Giapintzakis, D. 



Raman Scattering in Cuprate Superconductors 27 

M. Ginsberg and J.M. Mochel, Phys. Rev. B 44, 10158 B. Bucher et ai, Phys. Rev. 
Lett. 70, 2012 (1993); T. Ito, K. Takenaka, and S. Uchida, Phys. Rev. Lett. 70, 3995 
(1993); A. P. Mackenzie, S. R. Julian, C. T. Lin and D. C. Sinclair, Phys. Rev. B 53, 
5848 (1996). 

29. K. Kishio, in Coherence in Superconductors, edited by G. Deutscher and A. 
Revcolevschi (World Scientific, Singapore, 1995). 

30. D. N. Basov et al., Phys. Rev. B 50, 3511 (1994); C. C. Homes et al, Physica C 254, 
265 (1995); preprint. 

31. We note that we also are neglecting the important effects of doping and/or disorder on 
the pairing interaction. This undoubtedly plays a large part of the story concerning T c 
suppression, but does not affect the details of the overall lineshape. To proceed further 
requires a microscopic mechanism for pairing itself. 

32. P. J. Hirschfeld, J. Phys. Chem. Solids 56, 1605 (1995). 

33. A. A. Abrikosov and L. P. Gorkov, Pis'ma Zh. Eksp. Teor. Fiz. 40, 351 (1994) [Sov. 
Phys. JETP 40, 1155 (1985)]. 

34. P. W. Andcrson, J. Phys. Chem. Solids 11, 26 (1959). 

35. T. P. Devereaux and A. P. Kampf, unpublished. 

36. B. Stojkovic and D. Pines, Phys. Rev. Lett. 76, 811 (1996); to appear in Phys. Rev. B. 

37. A. V. Chubukov, D. K. Morr, and K. A. Shakhnovich, Phil. Mag. B 74, 563 (1996); D. 
Pines, preprint. 



